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HIGHLIGHTS 


•  Local  current  density  is  measured  by  two  specially  designed  single  cells. 

•  Variations  of  local  current  density  during  charge/discharge  process  are  presented. 

•  Effect  of  electrolyte  flow  rate  on  current  density  variation  is  shown. 

•  Lower  limit  of  state  of  charge  during  discharge  process  is  evaluated. 
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This  article  presents  a  preliminary  study  of  the  measurement  of  local  current  density  in  all-vanadium 
redox  flow  batteries.  Two  batteries  are  designed  and  manufactured  in  this  study,  and  the  experi¬ 
mental  results  are  compared.  In  the  first  cell,  the  current  collector  is  divided  into  25  segments,  and  the 
flow  field  plate  is  not  segmented,  whereas  in  the  other  cell,  the  flow  field  plate  is  segmented.  The  effects 
of  the  electrolyte  flow  rate  on  the  battery  efficiencies  and  the  local  current  density  variation  are 
investigated.  The  experimental  results  show  that  the  current  density  near  the  outlet  significantly  de¬ 
creases  when  the  discharge  capacity  approaches  zero.  In  addition,  the  battery  has  a  larger  discharge 
depth  at  a  higher  electrolyte  flow  rate. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 


At  the  negative  electrode: 


All-vanadium  redox  flow  batteries  (VRFBs)  have  garnered  sub¬ 
stantial  attention  over  the  past  few  years  because  of  renewable 
energy  storage  requirements.  A  VRFB  discharges  and  stores  elec¬ 
tricity  through  the  redox  reaction  between  vanadium  ion  oxidation 
states.  In  a  VRFB  system,  V2+/V3+  and  V4+/V5+  redox  couples  serve 
as  the  negative  and  positive  electrolytes,  respectively  1,2  .  During 
discharge,  V2+  is  oxidized  to  V3+,  while  V5+  is  reduced  to  V4+.  In  the 
charging  process,  these  reactions  proceed  in  the  opposite  direction. 

At  the  positive  electrode: 

Discharge  ~ 

VOj  +  2H+  +  e~  <  *  V02+  +  H20.  (1) 

Charge 
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Discharge 

Charge 


V3+  +e“. 


Both  the  negative  and  positive  electrolytes  are  stored  in  reser¬ 
voirs  outside  the  VRFB  cell  and  are  pumped  into  the  cell  stack 
during  electrochemical  reactions.  Because  of  this,  the  energy  ca¬ 
pacity  can  be  increased  by  simply  adding  more  electrolytes  to  the 
reservoirs,  and  the  output  power  is  determined  by  the  cell  number 
and  active  area  of  the  electrodes.  Because  the  energy  capacity  is 
determined  by  the  amount  of  electrolyte  stored  in  the  reservoirs 
and  because  the  power  output  is  determined  by  the  cell  stack,  the 
power  and  energy  of  a  VRFB  system  can  be  individually  designed. 
Moreover,  the  positive  and  negative  electrolytes  are  separately 
stored,  resulting  in  low  self-discharging.  Because  of  these  advan¬ 
tages  over  other  energy  storage  systems,  VRFB  systems  are 
considered  potential  candidates  for  large-scale  energy  storage  and 
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Power  Supply 


VRFB  Cell 


Electrolyte  Tank 


Fig.  1.  A  schematic  of  the  experimental  setup. 


have  been  coupled  to  intermittent  renewable  power  sources,  such 
as  solar  or  wind  power  systems  [3-6].  In  the  coupled  system,  VRFB 
systems  perform  load  leveling  and  peak  shaving. 

The  major  challenges  for  the  commercialization  of  VRFBs  are  the 
key  materials  and  system  optimization,  which  significantly  influ¬ 
ence  the  efficiency  of  a  VRFB  system.  Several  studies  have  investi¬ 
gated  the  key  materials,  including  electrolyte  additives  [7], 
sulfate-chloride  mixed  electrolytes  [8],  carbon  paper  electrode 
coated  with  supported  tungsten  trioxide  [9  ,  and  modified  sepa¬ 
rators  [10-17  ,  but  few  have  focused  on  the  system  efficiency 
[18-20]. 

Ma  et  al.  [18]  experimentally  studied  the  effect  of  the  elec¬ 
trolyte  flow  rate  on  the  performance  of  VRFBs  operating  at 
different  current  densities  and  proposed  an  optimal  strategy  for 
the  electrolyte  flow  rate.  They  suggested  that  at  each  operating 
current  density,  the  electrolytes  should  be  supplied  with  either 
low  or  high  flow  rates,  depending  on  the  cell  voltage.  Ma's 
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Fig.  3.  Charge-discharge  curves  of  cell  A  at  60  mA  cm 


experimental  results  also  showed  that  the  capacity  and  efficiency 
increased  with  increasing  electrolyte  flow  rate.  Benjamin  et  al. 
[19]  compared  the  battery  efficiencies  determined  using  cutoff 
voltages  and  SOC  limits.  Tang  et  al.  20]  modeled  the  concentra¬ 
tion  overpotential  as  a  function  of  the  electrolyte  flow  rate  to 
investigate  the  effect  of  the  electrolyte  flow  rate  on  battery  effi¬ 
ciency.  They  suggested  an  optimized  variable  flow  rate  with  a  flow 
factor  of  7.5,  which  resulted  in  higher  system  efficiencies  over 
constant  flow  rates.  When  the  electrolyte  flow  rate  increased 
because  of  a  low  reactant  concentration,  the  concentration  over¬ 
potential  was  reduced.  Although  the  effect  of  the  electrolyte  flow 
rate  on  cell  performance  has  been  studied,  its  effect  on  local 
current  density  has  not  been  reported  as  far  as  we  know. 

In  this  study,  two  VRFB  cells  are  designed  for  the  measurement 
of  local  current  density.  The  major  difference  in  the  design  is  the 
method  of  dividing  the  active  area.  The  experimental  results  of  the 
two  designs  are  compared  and  discussed.  The  effect  of  the  elec¬ 
trolyte  flow  rate  on  the  current  density  distribution  is  also  studied. 


2.  Experimental 

The  experimental  setup,  which  consists  of  the  VRFB  single  cell, 
electrolyte  recirculation  system,  and  the  measurement  system,  is 
illustrated  as  Fig.  1. 


outlet 


(b) 
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Fig.  2.  Schematics  and  photos  of  the  two  cell  designs  used  in  this  study,  (a)  Cell  A:  gold-plated  current  collector  divided  into  25  segments,  (b)  Cell  B:  graphite  flow  field  plate  divided 
into  25  segments. 
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Table  1 

Efficiencies  of  cell  A  at  60  mA  cur2. 


Operating  condition 

CE  (%) 

VE  (%) 

EE  (%) 

DC  (mAh  mU1) 

(A-i60-Q25) 

92.7 

76.0 

70.5 

17.34 

(A-i60-Q50) 

94.6 

79.0 

74.7 

24.28 

(A460-Q75) 

94.4 

79.3 

74.9 

26.42 

(A-i60-Q100) 

94.3 

80.0 

75.5 

27.31 

2.2.  VRFB  single  cell 


Fig.  5.  Charge-discharge  curves  of  cell  B  operated  at  40  mA  cm  2. 


In  this  study,  we  designed  and  manufactured  two  single  cells 
with  different  structures  for  the  local  current  density  measure¬ 
ment.  The  current  density  distribution  measurement  technique  for 
the  two  cells  was  similar.  The  active  area  of  each  cell  was  divided 
into  25  segments.  Each  segment  was  connected  with  a  50  mQ 
current  shunt. 

For  both  of  the  two  cells  with  an  active  area  of  100  cm2 
(10  x  10  cm),  a  proton  exchange  membrane  (Nation  117,  Dupont, 
USA)  was  employed  as  the  separator.  The  porous  electrodes 
(100  mm  x  100  mm  x  6  mm)  were  made  of  graphite  felt  (C0S1011, 
CeTech,  Taiwan).  The  major  difference  between  these  two  cells  is 
described  below. 

For  cell  A,  the  current  collector  of  the  negative  side  was 
replaced  by  25  gold-coated  copper  segments  (18  x  18  x  1  mm), 
as  shown  in  Fig.  2(a).  Each  segment  was  soldered  with  a  wire 
connected  with  a  current  shunt.  To  decrease  the  contact 


Table  2 

Efficiencies  of  cell  B  at  40  mA  cm  2. 


Operating  condition 

CE  (%) 

VE  (%) 

EE  (%) 

DC  (mAh  ml”1) 

(B-i40-Q50) 

95.3 

60.9 

58.1 

16.5 

(B-i40-Q100) 

94.5 

61.5 

58.1 

17.7 

resistance,  a  flexible  graphite  sheet  was  placed  between  the  flow 
field  plate  and  the  25  segments. 

For  cell  B,  a  space  (100  x  100  x  5  mm)  for  the  porous  electrode 
and  the  25  holes  (18  x  18  x  10  mm)  for  the  graphite  segments  were 
machined  on  an  end  plate  made  of  polyvinylchloride  (PVC),  as 
shown  in  Fig.  2(b).  Each  graphite  segment  was  plugged  with  a 
banana  connector.  Sealant  was  applied  between  the  graphite  seg¬ 
ments  and  the  end  plate  to  avoid  leakage  of  the  electrolyte. 


(a) 


Fig.  4.  Current  density  distributions  of  cell  A  at  (a)  the  beginning  of  the  charging  process;  (b)  the  end  the  charging  process;  (c)  the  beginning  of  the  discharging  process;  and  (d)  the 
end  of  the  discharging  process.  The  operating  condition  is  (A-i60-Q100). 
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Fig.  6.  Current  density  distributions  of  cell  B  at  (a)  the  beginning  of  the  charging  process;  (b)  the  end  the  charging  process;  (c)  the  beginning  of  the  discharging  process;  and  (d)  the 
end  of  the  discharging  process.  The  operating  condition  is  (B-/40-Q50). 


Both  the  inlets  of  the  anode  and  cathode  electrolytes  were  placed 
at  the  bottom  left  of  the  active  area,  and  the  outlets  were  placed  at 
the  top  right,  as  indicated  in  Fig.  2(a).  The  inlet  and  outlet  positions  of 
both  cells  were  the  same.  The  numbers  shown  in  Fig.  2(b)  help 
indicate  the  position  of  the  measured  local  current  density. 

2.2.  Electrolyte  system 

The  electrolyte  was  obtained  by  dissolving  1.5  M  VOSO4  in  2.0  M 
H2SO4.  The  theoretical  energy  capacity  is  approximately  40  mA  per 
ml  of  electrolyte.  Both  the  positive  and  negative  electrolyte  reser¬ 
voirs  were  initially  filled  with  the  V4+  solution.  The  positive  elec¬ 
trolyte  reservoir  was  initially  200  ml,  while  the  negative  side  was 
100  ml.  In  the  first  charging  step,  the  V4+  in  the  electrolyte  was 
converted  to  V5+  and  V2+  for  the  positive  and  negative  sides, 
respectively.  Next,  half  of  the  positive  electrolyte  was  removed  to 
make  the  electrolyte  volume  in  the  two  reservoirs  equal.  A 
continuous  nitrogen  flow  was  fed  to  the  negative  electrolyte 
reservoir  to  avoid  oxidation  of  the  V2+.  Two  diaphragm  pumps 
(SMART  digital  DDA7.5-16AR-PVC/V/C,  Grundfos,  Denmark)  were 
employed  to  recirculate  the  positive  and  negative  electrolytes  be¬ 
tween  the  reservoirs  and  the  VRFB  cell. 

2.3.  Measurement  system 

The  VRFB  cell  was  charged  with  a  power  supply  (U8001A, 
Agilent,  USA)  and  discharged  using  an  electronic  load 


(PLZ664WA,  Kikusui  Electronics,  Japan).  The  cell  voltage  and 
current  were  monitored  and  recorded  using  a  data  acquisition 
system  (NI  9205,  National  Instrument,  USA)  for  subsequent 
analysis. 

2.4.  Experimental  procedure 

The  current  density  was  recorded  when  the  VRFB  was 
operated  with  selected  electrolyte  flow  rates  and  load  currents. 
Both  the  charge  and  discharge  current  densities  were  the  same 
and  constant.  At  each  current  density,  both  the  positive  and 
negative  electrolytes  were  delivered  and  recirculated  with 
appropriate  flow  rates.  Under  each  operating  condition,  the  cell 
was  charged  and  discharged  between  0.7  and  1.7  V  for  three 
cycles,  and  the  experimental  data  of  the  third  cycle  was  used  to 
calculate  the  battery  efficiencies.  The  coulombic  efficiency  (CE), 
voltage  efficiency  (VE),  energy  efficiency  (EE),  and  discharge  ca¬ 
pacity  (DC)  were  calculated  using  the  following  equations, 
respectively: 

/  ^discharge 

CE  =  -  (3) 

J  ^'charge 

yg  =  ^avg,  discharge  ^ 

^avg, charge 
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Fig.  7.  Current  density  distributions  of  cell  B  at  (a)  the  beginning  of  the  charging  process;  (b)  the  end  the  charging  process;  (c)  the  beginning  of  the  discharging  process;  and  (d)  the 
end  of  the  discharging  process.  The  operating  condition  is  (B-/40-Q.100). 


/  ^VRB, discharge 

EE  =  J— -  (5) 

/  ^VRB, charge 

DC  =  J  i'discharge(0^-  (6) 


3.  Results  and  discussion 

Although  the  experimental  results  indicate  that  cell  A  is  not 
suitable  for  the  measurement  of  local  current  density,  some  data 
from  cell  A  under  selected  operating  conditions  are  presented  and 
discussed.  The  operating  conditions  are  described  using  the  nota¬ 
tion  (A/B-i#-Q#),  where  A/B  denotes  cell  A  or  B,  i  represents  the 
current  density  in  mA  cnrT2,  and  Q.  is  the  electrolyte  flow  rate  in 
ml  min-1. 

3.1.  Current  density  distribution  of  cell  A 

Fig.  3  shows  the  charge-discharge  curves  of  cell  A  while  oper¬ 
ating  at  60  mA  cm-2  with  four  electrolyte  flow  rates:  25, 50,  75,  and 
100  ml  min-1.  The  battery  efficiencies  of  cell  A  are  compared  in 
Table  1  based  on  the  experimental  data  in  Fig.  3.  This  shows  that  the 


DC  increases  with  increasing  electrolyte  flow  rate.  This  tendency  is 
similar  to  an  ordinary  VRFB  18].  Four  selected  snapshots  of  the 
current  density  distribution,  corresponding  to  the  marks  in  Fig.  3, 
are  presented  in  Fig.  4. 

Fig.  4(a)-(d)  represents  the  current  density  distributions  during 
the  initial  charging  period,  the  final  charging  period,  the  initial 
discharging  period,  and  the  final  discharging  period,  respectively.  It 
can  be  seen  that  the  current  density  distribution  does  not  vary 
significantly  with  time.  This  may  be  caused  by  the  contact  re¬ 
sistances  among  the  25  segments  not  being  equal,  resulting  from 
non-uniform  compression  pressure  on  each  segment.  As  a  result, 
the  current  density  distribution  is  mainly  affected  by  the  contact 
resistance.  Another  possible  reason  is  the  in-plane  current  flow  in 
the  5  mm-thick  flow  field  plate,  resulting  in  a  redistribution  of 
current  density  when  current  flows  from  the  electrode  side  to  the 
current-collector  side  through  the  flow  field  plate. 

3.2.  Current  density  distribution  of  cell  B 

Cell  B  is  designed  to  eliminate  the  effects  of  the  in-plane  current 
flow  on  the  current  density  distribution.  The  charge-discharge 
curves  of  cell  B  operated  at  40  mA  cm-2  are  shown  in  Fig.  5,  and  the 
battery  efficiencies  are  listed  in  Table  2.  Compared  to  cell  A,  cell  B 
shows  significantly  lower  VE  and  DC.  This  may  be  because  the 
contact  resistance  between  the  banana  plug  and  graphite  is  higher 
than  that  between  the  current  collectors  and  the  graphite  plate  in 
cell  A. 
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Fig.  8.  Variations  of  the  local  current  density  with  respect  to  the  SOC  of  cell  B  operated 
at  (B-/40-Q50)  during  the  (a)  charging  and  (b)  discharging  process. 

Figs.  6  and  7  show  the  current  density  distributions  of  cell  B 
operated  under  the  conditions  (B-i40-Q50)  and  (B-i40-Q100), 
respectively.  Fig.  6(a)  and  (b)  shows  that  during  the  charging 
process,  no  considerable  variation  in  current  density  distribution 
is  observed.  However,  at  the  end  of  discharging  process,  the 
current  density  close  to  the  outlet  decreases  significantly  as 
shown  in  Fig.  6(c)  and  (d).  Because  the  concentration  of  V2+  (and 
V5+)  in  the  electrolyte  is  the  lowest  at  the  outlet,  the  local  per¬ 
formance  of  the  cell  close  to  the  outlet  is  the  lowest.  It  can  also 
be  seen  in  Fig.  6(d)  that  at  the  end  of  discharging  process,  the 
local  current  density  of  the  right  half  area  of  the  cell  is  smaller 
than  that  of  the  left  half.  This  may  be  attributed  to  the  electrolyte 
flow  distribution.  In  cell  B,  the  inlet  is  placed  lower  than  the 
bottom  of  the  active  area  as  depicted  in  Fig.  2(b).  Most  of  the 
inflow  electrolyte  initially  goes  upward  and  then  toward  the 
right,  resulting  in  a  higher  concentration  of  V2+  (and  V5+)  in  the 
left  half  of  the  cell.  The  electrolyte  flow  distribution  can  be 
improved  through  the  flow  field  design  of  the  frame  [21  ;  how¬ 
ever,  its  effect  on  current  density  distribution  is  not  focused  at 
current  stage  of  this  study. 
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Fig.  9.  Variations  of  the  local  current  density  with  respect  to  the  SOC  of  cell  B  operated 
at  (B-i40-Q100)  during  the  (a)  charging  and  (b)  discharging  process. 


Fig.  8  shows  the  variations  of  the  local  current  density  with 
respect  to  the  DC  of  cell  B  operated  under  (B-i40-Q50)  during  the 
charging  and  discharging  processes.  During  the  charging  process, 
the  local  current  density  exhibits  a  stable  distribution;  however, 
during  the  discharging  process,  a  variation  of  the  current  density 
distribution  occurs  when  the  DC  approaches  approximately  20%. 
Furthermore,  the  local  current  density  of  a  number  of  segments 
close  to  the  outlet  shows  dramatic  decreases  at  a  DC  of  10%.  Some  of 
them  even  approach  zero.  In  contrast,  the  local  current  density  of  a 
number  of  segments  close  to  the  inlet  increases  during  a  galvanic 
discharge.  It  has  been  experimentally  shown  that  the  battery  effi¬ 
ciency  decreases  with  increasing  current  density.  Thus,  the  lower 
limit  of  the  DC  during  the  discharging  process  can  be  determined  as 
10%  under  the  operating  condition  (B-i40-Q50). 

The  variations  of  the  local  current  density  at  an  electrolyte  flow 
rate  of  100  ml  min1  are  shown  in  Fig.  9.  Fig.  9(b)  suggests  that  the 
lower  limit  of  the  DC  is  approximately  5%.  Equivalently,  the 
discharge  depth  is  increased  by  increasing  the  electrolyte  flow  rate. 
However,  in  contrast,  higher  electrolyte  flow  rates  result  in  larger 
power  consumptions  by  the  pumps.  Therefore,  the  effect  of  the 
electrolyte  flow  rate  on  the  local  current  density  distribution  and 
the  discharge  depth  must  be  evaluated  to  determine  an  optimal 
electrolyte  flow  rate.  In  the  future,  the  effects  of  the  electrode 
design  on  the  local  current  density  will  be  investigated  using  a 
newly  designed  cell. 

4.  Conclusion 

To  measure  the  current  density  distribution,  two  VRFB  cells  with 
an  active  are  of  100  cm2  were  designed  in  this  study.  In  the  first  cell, 
the  gold-plated  current  collector  of  the  negative  side  is  divided  into 
25  segments.  In  the  other  cell,  the  positive  and  negative  graphite 
flow  field  plates  are  divided  into  25  segments.  It  is  shown  that  the 
cell  design  significantly  influences  the  experimental  results.  The 
accuracy  of  the  measurement  may  be  attributed  to  both  the 
thickness  of  the  flow  field  plates  and  the  contact  resistance.  The 
second  cell  with  divided  graphite  field  plates  exhibits  a  significant 
improvement  in  the  measurement  of  the  local  current  density. 
However,  its  performance  decreases  because  of  the  increased 
contact  resistance.  The  experimental  results  show  that  the  current 
density  close  to  the  outlet  significantly  decreases  when  the  DC 
approaches  zero.  Moreover,  the  discharge  depth  is  larger  at  a  higher 
electrolyte  flow  rate. 
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